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Abstract 

We analyze the indirect astrophysical signatures of secluded models of WIMP dark mat- 
ter, characterized by a weak-scale rate for annihilation into light MeV-scale mediators which 
are metastable to decay into Standard Model states. Such scenarios allow a significant 
enhancement of the annihilation cross section in the galactic halo relative to its value at 
freeze-out, particularly when the mediator is light enough for this process to proceed through 
radiative capture to a metastable WIMP-onium' bound state. For MeV-scale vector me- 
diators charged under a hidden U(l)' gauge group, the enhanced annihilation rate leads 
predominantly to a sizable excess positron flux, even in the absence of astrophysical boost 
factors. 
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1. Introduction 



The growing level of astrophysical and cosmological evidence for dark matter has led in re- 
cent years to an expansion in experimental programs that aim to detect its non-gravitational 
interactions pQ. These range from direct production at colliders to the recoil of galactic dark 
matter on nuclei in underground detectors, and indirect detection of annihilation products 
from various parts of the galactic halo. The primary theoretical motivation for anticipating 
interactions with the Standard Model (SM) at a detectable level is the apparent simplicity 
of the thermal WIMP paradigm. Namely, the observation that a thermal relic with weak 
scale mass and annihilation cross-section naturally provides roughly the correct cosmological 
abundance [2]. However, even this additional theoretical input still leaves a vast range of 
possibilities, and thus it is far from clear which of the various experimental or observational 
strategies will have the best sensitivity. It is therefore important to ensure that the obser- 
vational reach is as broad as possible, given reasonable theoretical expectations, and to that 
end it is useful to explore the differing levels of sensitivity for various detection strategies 
within the simplest consistent models for thermal dark matter. 

A generic problem is that while the thermal WIMP scenario imposes constraints on the 
annihilation cross-section, it clearly does not directly constrain the scattering cross-section 
with matter, and it is possible to consider scenarios where these rates differ significantly. 
An example is the class of secluded WIMPs j3], or alternatively a dark hidden sector, where 
WIMPs annihilate with a weak-scale rate to metastable mediators which are in turn very 
weakly coupled to the Standard Model. In such scenarios, the direct scattering cross-section 
is negligible, and the primary detection strategy is via observation of the products of anni- 
hilation in parts of the galactic halo. Various annihilation products may be observed, i.e. 
7-rays, anti-matter (positrons and anti-protons), and also neutrinos from annihilation in the 
sun and elsewhere, depending on the form of the model. Thus, it is important to know which 
if any of these annihilation products would provide the best sensitivity given current and 
future observational capabilities. Such an investigation is timely given that existing data 
from the HEAT [I] and AMS-01 [5] experiments indicating a possible excess in the relative 
positron fraction n e +/(n e - + n e +), and thus a possible enhancement in the production of 
O(10 GeV) positrons in the galactic vicinity of the solar system, recently received support 
from preliminary results from PAMELA [B] . Although such an excess - even if confirmed to 
be above the expected background level produced via cosmic rays [7] - could be due to vari- 
ous astrophysical sources, e.g. pulsars in the local neighbourhood, the recent announcement 
has generated a flurry of theoretical activity addressing the question of whether a positron 
excess could be tied to WIMP annihilation [8|. This recent work has built upon several ear- 
lier analyses [HJ [TO] , and a fairly general conclusion is that an excess at the level apparently 
observed by HEAT and AMS-01 (and PAMELA) would require a nontrivial enhancement of 
the annihilation rate in the galaxy relative to its rate at freeze-out for a thermally populated 
WIMP. Such enhancements are often conjectured (in the form of boost factors) to result 
from local features of the halo distribution, but more generally this raises a challenge to see 
whether generic WIMP scenarios could produce an anti-matter flux at a level sufficiently far 
above background to be observed. 

In this note, in the context of secluded WIMP models, we point out that if the mass of 
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the mediator m v satisfies the following bounds, 



2m e <m v < p DM ^ m DM a' 



(1) 



where pdm = wdm^ is the typical nonrelativistic WIMP momentum inside the galactic halo, 
and a' is the strength of the interaction in the dark sector, the resulting annihilation rate is 
enhanced over the cosmological rate by a factor of 0(a'/v) > 1, while the decay products 
consist mostly of electron and muon pairs and photons. Most importantly, if my is less than 
the 'WIMP-onium' binding energy -Ebind, i-e. if 



then annihilation inside the halo proceeds via the intermediate formation of a metastable 
WIMP-onium bound state, and the enhancement factor over the cosmological annihilation 
rate grows to O(10 a' /vdm)- Interestingly enough, this appears to be of sufficient magnitude 
to boost the observed positron fraction in the multi-GeV range to a level that would be 
observed as an excess above background without the need for astrophysical boost factors. 

To provide a concrete example, we calculate the indirect detection signal in a model of 
secluded dark matter with a vector U(l)' mediator V kinetically mixed with the photon, 
noting that similar results can be obtained for scalar mediation. In the next section, we 
discuss the various annihilation processes, and the distinctions between those most relevant 
for co mo logical freeze-out, and for producing observable signals in the galaxy today. In 
Section 3, we briefly discuss the observational implications of galactic annihilation focusing 
on excess positrons and 7's, and conclude in Section 4 with some other implications of this 
dark matter scenario. 



The characteristic feature of secluded WIMPs [3] is that the WIMP in the hidden sector 
communicates with the Standard Model via a mediator, with mwiMP > ^mediator, such that 
direct annihilation into a pair of mediators is always possible. The correct relic density may 
then be achieved while imposing only very mild constraints on the coupling of the mediator 
to the SM - a mediator lifetime of up to a second, so that decays occur before primordial 
nucleosynthesis (BBN), is sufficient - and this allows the WIMP to be secluded from the SM 
with a very small cross-section for direct scattering. This generic feature leads us to focus 
on the annilation process as a means for indirect detection of such dark matter scenarios. 

The annihilation details are sensitive to the relevant mass scales of the WIMP and the 
mediator, and scenarios in which the mediator may couple to the SM at the renormalizable 
level are of particular interest since the mediator can then naturally be rather light. There 
are only a limited number of options for introducing such couplings, often referred to as 
portals, and most of them have been considered at length in a number of publications 
[TTJ [12]. One of the most natural possibilities that leaves the mass-scale of the mediator 
as a free parameter is the kinetic mixing portal between the hypercharge gauge boson 
and a U(l)' mediator [13]. The only renormalizable coupling to the SM occurs through 
kinetic mixing, £ mix ~ kV^B^, and the SM is neutral under U(l)'. This is the scenario we 
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2. Annihilation of Secluded WIMPs 
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Figure 1: WIMP annihilation: on the left, the direct decay to two metastable on-shell Vs, which ultimately 
decay to SM states; and on the right, the decay through the formation of a metastable (para) WIMP-onium 
state, which occurs with an enhanced rate for non-relativistic WIMPs in the galactic halo. 



will explore here, as the vector nature of the mediator has important consequences for the 
relative fluxes of different annihilation products. 

We assume that the hidden sector is a singlet under the SM gauge group, while carying a 
charge under U(l)'. We further assume that U'(l) is spontaneously broken, and restricting 
the discussion to renormalizable couplings, choose a Lagrangian of the following form: 

X 1 2 2 

^WIMP+V-portal = AviMP _ l^V^li + ^v^^F^ + Cy- (3) 

Here we have retained the mixing of V with the photon field strength F^, and combined the 
residual Higgs' terms in Ch>- The simplest form of the WIMP Lagrangian is either a Dirac 
fermion or complex scalar charged under U(l)', 

^wimp = i>{ iD iSYn ~ m ^ + C(Am^) fermionic DM, (4) 
£wimp = (D^)*(D^) - ml4>*4> + C(A m4) ) bosonic DM, (5) 

where = <9 M + ie'V^ is the usual covariant derivative in terms of the U(l)' gauge coupling 
e'. The WIMPs ip or <fi will be secluded provided m v < m^u), and for the present discus- 
sion the relevant regime will be my <C m^uy In general, one can introduce mass terms 
£(Am^j) that lift the degeneracy between ip, and their charge-conjugated copies. To be 
consistent with U(l)' gauge invariance the mass splitting has to be proportional to the scale 
of spontaneus symmetry breaking, Am^ ~ v' — mv/e', which in this paper we choose to be 
small. While we note that such a splitting may introduce additional interesting signatures 
for direct and indirect detection [TJJ [151 ES]> we leave this issue aside here and treat ip and 
(f) respectively as a pure Dirac fermion and a charged scalar. 



2.1 Relic Abundance 



The dominant annihilation process for determining the relic abundance was discussed in [3]. 
Provided my < m^, the two-step process takes the form (see Fig. [1]): 

(1) V + i> -> on-shell V + V, 

(2) V -> virtual 7, Z -> SM states. 
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The cosmological abundance is determined only by process (1), and the spin-averaged cross 
section for %p + %l)~ annihilation into pairs of (unstable) V bosons is then given by 



.,_ n(a') 2 I ml m v ^m^ n(a') 2 

while for <p + (j)~ annihilation the result of ([6]) should be multiplied by 2. In order to use (EI) 
in computing the relic cosmological abundance, one has to average it over the three different 

collision types: H — , ++ and , where + and — refer to the U(l)' charge, so that in 

the absence of a significant primordial asymmetry between particles and antiparticles in the 
WIMP sector, (av) = \ av \ 2 y ■ Assuming that ip (or 0) is the dominant component of dark 
matter, and relating the cross section to the measured cold dark matter energy density, we 
obtain: 

i < Q DM h 2 « 0.1 => (av) = 2.5 x 10- 26 cm 3 s~\ (7) 

V 9*(Tf) x (av) 

very close to the conventional order-of-magnitude estimate. Here Xf is the inverse freeze-out 
temperature in units of the WIMP mass (xf = m^/Tf ~ 20), and g* ~ 0(100) is the number 
of effective degrees of freedom. Eq. (j7J implies that in the limit of small mixing, «<1, the 
correct dark matter abundance is achieved if the WIMP masses and couplings satisfy the 
relation(s) 

a'^10- 2 x( _ m t A o^lO" 2 *' '"• 



270 GeV/ ' V380 GeV 

(These formulae also correct a small numerical error in Eq. (13) of Ref. [3].) The conditions 
([8]) are easily satisfied for a rather natural range of m^u) and a'. The constraints that arise 
on the mixing parameter k, which determines the ultimate decay rate of V — > SM states, 
are very mild and k can be taken very small indeed. As discussed in [3], the only constraints 
one has to impose are that the decay of V (and also h') occur before the start of Big Bang 
nucleosynthesis (BBN). Given m' h > my/2, only the decays of V are sensitive to the mixing: 

r >'^ ^ " 2 (low) < 9 > 

If the mediators are light, then my is also constrained by BBN even if the V bosons remain 
in thermal equilbrium. An analysis of BBN in the presence of unstable MeV-scale relics 
coupled to the electron-positron plasma [17] shows that relics with mass above 4 MeV are 
consistent with observations, and thus we adopt this value as a lower bound in our analy- 
sis, my ^ 4 MeV. Requiring in addition that V decays remain in thermal equilibrium, i.e. 
Ty > Hubble Rate [T ~ 0.05m,/,] at freeze-out, would ensure the initial thermal and chemical 
equilibirium for WIMPs used in the derivation of the abundance formula, leading to a signif- 
icantly tighter constraint on k [3J . However, this may be relaxed back to (Q if some new UV 
physics ensures proper thermal contact between the dark matter and SM sectors at higher 
temperatures. For the rest of this paper we assume that k, e'/e, so that its precise value 
will not play a role in calculations of the WIMP abundance and astrophysical signatures of 
WIMP annihilation. To define a fidicual parameter range, we shall mostly consider coupling 
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constants a' on the order of the SM gauge couplings at the weak scale, and take the dark 
matter mass to lie in a rather liberal range, 100 GeV < m^u,) < few TeV. 

The recoil signal in direct detection depends rather sensitively on all the parameters of 
the model: a', my, k and m^). The relevant WIMP-nucleus scattering cross sections have 
been calculated previously in [31 HE], which for my > [iv - where \i is the reduced mass and 
v the relative velocity of the WIMP-nucleus system - is mediated by the electromagnetic 
charge radius of the WIMP: r 2 = 6ne' /(em v ). The cross section per nucleon in this limit is 

(Z\ 2 16tt K 2 aa'm 2 v , 

'=U) X — ST-*- (10) 

For my of electroweak scale, small values of the mixing parameter k strongly suppress the 
scattering cross section relative to its 'natural' weak-scale value, implying that the WIMP 
is indeed 'secluded'. However, the limit of small my considered in this paper may entail 
significant direct-detection sensitivity to the parameter space of the model owing to the my 4 
scaling of (1101) . The current experimental sensitivity to the WIMP-nucleon cross section, 
a ~ 10~ 43 cm 2 , in the most favorable range of w^u) ~ 100 GeV, probes the following 
combination of masses and couplings, 

2 f a' \ /100 MeV\ 4 17 

K x (hfO(^H ~ fewxl ° ■ 

which is quite close to the lower bound on k in Eq. (Q. Note that for my < fiv this formula 
needs to be modified to correct for the dependence on the energy of the scattering particles. 
Although (ITTj) suggests impressive sensitivity to k, it should be emphasized that even a 
modest splitting Am^^) ~ O(MeV) in the dark sector removes this completely as r c is zero 
for the Majorana fermions or real scalars [IB] , while off-diagonal transitions are kinematically 
forbidden for v ~ 10~ 3 . As a consequence, much larger values of k can still be accomodated 
for small my. 

2.2 Enhanced annihilation inside the halo 

Before we present the annihilation rates in secluded WIMP models, we would like to briefly 
comment on the generic mechanisms for enhancing the annihilation rate in the halo relative 
to the freeze-out value of ~ 3 x 10 _26 cm 3 s _1 . We introduce a dimensionless parameter, 

Af = ( 7 )hal ° , (12) 

freeze— out 

that in the Born approximation is expected to be of order one if the annihilations are pre- 
dominantly s-wave, and much less than one for p-wave annihilation. This expectation is 
based on the fact that the halo velocities of WIMPs are much smaller than the correspond- 
ing velocities at the early cosmological freeze-out epoch. For the purposes of this paper it is 
sufficient to assume an approximately Maxwellian velocity distribution inside the halo with 
'temperature' T h and introduce a parameter Xh, 

ii, = ?~3x 10 6 , (13) 
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that corresponds to the choice of an r.m.s. WIMP velocity of Vh ~ 200 km/s. Below we list 
the important exceptions that allow for M 1: 

• A: Long-range attractive interaction between WIMPs. In this case, the enhancement 
is given by the familiar Sommerfeld factor that accounts for the modification of the 
scattering wavefunction, relative to the leading-order Born approximation, in the pres- 
ence of an attractive force. We will assume: (i) that the annihilation is s-wave; (ii) 
an exaggerated hierarchy of scales in Eq. ([T]), namely m v l 3> ob and my <C Ap^ 
where as ~ (c/m^,/2) _1 is the 'Bohr radius' and Ap^ ~ m^Vh is the relative WIMP 
momentum. This then implies as a consequence that a'/vt ^> 1, and the enhancement 
factor is 



For simplicity we have disregarded any modifications of the freeze-out rate, assuming 
that a'/vf <C 1. 

• B: Accidental near-threshold resonances. Very narrow, near-threshold, resonances can 
enhance the halo annihilation rate quite significantly. Such resonances may be due 
to the accidental presence in the particle spectrum of a mediator with the quantum 
numbers of the WIMP pair and a mass my = 2m^. It is more natural, however, to 
expect a resonance in the form of a bound state of two WIMPs, which is possible if 
there is an attractive force acting in the WIMP sector. It is not difficult to estimate 
the maximum possible enhancement in this case. It is well-known that if the [ther- 
mal] spread in energy of interacting non-relativistic particles is much larger than the 
resonance width, the corresponding cross section can be written as a delta-function in 



where r res is a combination of the total, incoming, and outgoing widths as well as a 
multiplicity factor g, T res = gTi n r ou t/r to t, and E is the WIMP energy in the c.o.m. 
frame. If the resonant energy is within the energy distribution, the ratio of the two 
thermal rates is maximized and can be as large as 



Should all the requisite conditions be satisfied, a near-threshold resonance may indeed 
result in such a dramatic enhancement although in general the impact on the freeze- 
out rate, that we have ignored here, may not be negligible |19J . However, from the 
model-building perspective the probability of a narrow resonance right at the E = 
threshold appears rather low. Indeed, for m^ in the TeV range the maximum resonant 
enhancement occurs when the resonant level is within ~ 1 MeV of the threshold, which 
would be a rather special point in the {m^, m v , a'} parameter space. To have a bound 
state in the atractive Yukawa potential, one would typically expect the interaction 
range m v to be on the order of or larger than (a'm^)" 1 , and the heavier my and 
are, the more severe the fine-tuning will become. 




(14) 



energy, 




(15) 




(16) 
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C: New annihilation channels. While the two previous mechanisms have already been 
discussed at some length in the literature [2U1 EE], a third possibility - namely the 
enhancement due to new annihilation channels - has not been widely noticedfj] Specif- 
ically, we refer to the new recombination process (see Fig. 1), 

recombination : if) + ip — > (-0 — onium) + V, (17) 

(and similarly for <ft) which is kinematically open even in the limit — > if the 
condition (T5]) is satisfied. The subsequent fate of the ?/>-onium state is very different 
within the early Universe during freeze-out as compared to the galactic environment. 
In the halo, every ^-onium that is formed via the process (j!7p decays further to two or 
three ^/-bosons. During freeze-out, however, the annihilation rate of ^-onium into V's 
is strongly inhibited by thermal break-up, ^-onium + V — > 2ip. One can easily show 
that in the latter case the efficiency of annihilation, Br = r annih /(r annih + r broak _ up ), is 
much smaller than one. Thus, effectively only when the temperature drops below the 
binding energy does the process ( fTTl) serve as a new annihilation channel and, as we 
are going to see shortly, indeed dominate the annihilation rate in the galactic halo. 



When the recombination process ffT7j) is kinematically allowed according to (T5]) - which 
for m^ ~ 500 GeV and a' fixed according to ([8]) requires my < 50 MeV - the rate can be 
computed by generalizing the corresponding calculation for positronium to a finite vector 
mass my. Retaining only the direct recombination to the ground state of WIMP-onium, we 
arrive at a recombination rate that is independent of the spin of the annihilating particles. 
Assuming once again that my <C Ap^^ ~ m^^Vh, which is the range of parameters relevant 
here, we find 



3exp(4)mJ w V 2 J \ v 

where vy = [1 — Amy / (a'm^)) 2 ] 1 ! 2 is the velocity of the emitted V boson, which we assume 
is not significantly different from 1. Similar WIMP recombination processes were considered 
previously is Refs. [T6J [23]. The ^(0)-onium state will eventually decay to either 2V or 3V, 
with branching ratios of 25 and 75% for the fermionic WIMPs, and to 2V for bosonic dark 
matter. Since is large, the decay rates are rather fast, e.g. T^y ~ (a') 5 m^, followed by 

the subsequent much slower decays of V — > SM states with lifetimes controlled by the small 
mixing parameter k 2 . 

The crucial observation is that this recombinant annihilation process is significantly en- 
hanced in the galactic halo relative to the direct decay to 2V which dominates at freeze-out. 
Besides the expected a'/v factor, there is an additional well-known enhancement by the large 
numerical coefficient in ([18]) . Comparing (fl8l) to ([6]), we find the enhancement factors for 
the cases of fermionic and bosonic dark matter: 

A/^2o/-\ , A^~lo/-\ . (19) 



v / h \ v i h 



1 Recombination through a Coulomb-like force for a sub-dominant component of dark matter was consid- 
ered in [HI. 
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For a < a' < a w this implies a large enhancement factor of (9(100), as (v ~ l )h ~ (2x/ l /7r) 1 / 2 ~ 
1.5 x 10~ 3 . Moreover, this enhancement will occur as long as the conditions (OQ) and ([2]) are 
satisfied, and this does not require any tuning of resonant levels. The inclusion of recombina- 
tion to 2s WIMP-onium states, which is allowed if 16my < (a') 2 m^, would further increase 
this number by ~ 20%. If the emission of V bosons in recombination is not kinematically 
allowed, i.e. if my > -Ebind; but the conditions in ([1]) are still satisfied, the enhancement of 
e.g. the ip + ip — > 2V rate is given simply by Eq. (TH|) . which for the same value of a' is a 
factor of ~7(3) times smaller than J\f^^ in f[T§l) for fermionic(bosonic) WIMPs. 

The final step in the decay chain is the decay of the V bosons resulting from WIMP 
annihilation (possibly within a WIMP-onium bound state). Depending on the value of my, 
one may have different SM final states. In a generic secluded regime, my < m^h), such 
final states may involve leptons, both light and heavy quarks, and W pairs. Consequently, 
one should expect a non-negligible fraction of anti-protons being created in the process of 
hadronization of the decay products. However, in the domain of parameters where the halo 
annihilation rates are most significantly enhanced, my < -Ebind, only decays to light charged 
mesons and electron and muon pairs are kinematically possible, as -Ebind ^ 1 GeV even for 
multi-TeV WIMPs. As a consequence, all SM decay products will be highly boosted, and 
thus can provide attractive targets for indirect detection as 7-rays and positrons with energies 
of (9(100) GeV have a rather limited scope for production through conventional astrophysical 
sources. The final parameters of interest in this scenario are the relative branching fractions 
for the decay of V, which occurs through mixing with off-shell photons. As it turns out, the 
vector nature of the mediator in this case leads to a dominant leptonic branching fraction 
(e.g. V —>■ 7* — > e + e~), as direct decays to 27 are forbidden and the lowest order photonic 
decay rate is V — > 37 proceeding via a loop diagram. Direct decays of V's to neutrinos are 
also highly inhibited. Since all light charged mesons and muons eventually cascade down 
to electrons, positrons and neutrinos, we conclude that on average between two-to-three 
positrons are produced for each WIMP annihilation, while the yield of 7-quanta, resulting 
from internal bremsstrahlung in ^/-decays, is expected to be suppressed by an 0(a/ir) factoifl 
The dominant leptonic branching fraction has interesting observational consequences that 
we will explore in the next section. 

3. Observational Signatures 

The dominant observational signals for DM models of this type comprise an excess of high- 
energy 7-rays and positrons from annihilation. Since astrophysical sources generally fall as 
power-laws, the higher energy end of the spectrum, E ~ few 100 GeV, is the region with the 
best signal/background ratio. It also means that for the class of secluded models considered 
here, the 7-ray production rate (at least locally) will be of less importance than the positron 
rate, because of the ~ a/n suppression mentioned in the previous subsection. 

Turning to the observable flux, it is worth recalling that the signatures for both 7- 

2 We note in passing that for secluded WIMP models with scalar mediators in the O(100 MeV) range, the 
branching fraction of 7 quanta in the decay products will be significantly enhanced and can indeed dominate 
over other decay channels. 
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rays and positrons are highly dependent on details of the galactic halo, but in different 
ways. The observed 7-ffux is proportional to the line-of-sight integral, Jj s dl(a~ f v)t l p(l) 2 , 
where p(l) is the dark matter density along the line of sight. The flux thus samples the 
entire halo distribution, and is enhanced in directions where the density is largest, e.g. the 
galactic center. It is therefore sensitive to the fact that the inner halo profile is still not 
well constrained by simulations. Recently, several high energy sources believed to be of 
astrophysical origin have been detected in the direction of the galactic center by HESS [21] 
and EGRET [23], which unfortunately enhances the background in directions where the 
signal should be at its peak. In contrast, the observed positron flux, at least for energies 
near the peak of the decay spectrum, is a more 'local' signal. The observed flux is affected 
significantly by energy degradation through various loss mechanisms, such as synchrotron 
radiation and inverse Compton scattering. The flux can be represented in the form, 



where p$ is the constant normalization of the DM density in the local region, and dN e + / dE is 
the annihilation spectrum. For the 2V final state, the emitted V bosons are monoenergetic, 
which creates a flat spectrum of e + e~ pairs in the 'halo restframe', dN e +/dE = N Q(m^ — E). 
The propagation function Fp depends on the diffusion length Xd(E, E- m ), which determines 
the energy degradation from the injection energy E\ n to the observed energy E, and in 
general many other details of the halo distribution (see e.g. [26]). However, due to this 
energy degradation, if we focus only on the high-energy part of the spectrum where E ~ 
this effectively samples just the local neighbourhood where the dark matter density should 
be roughly constant. For E ~ ~ 500 GeV, the diffusion distance is probably at most 
1 kpc. In this approximation, F is obtained by integrating just over a local sphere and 
thus, while a full solution of the diffusion equation [27] will introduce some sensitivity to 
the halo profile, the high-energy part of the positron flux is primarily sensitive only to the 
local DM number density, as is also the case for direct-detection. Although p is itself quite 
uncertain, e.g. due to halo substructure, this effectively limits the DM-related uncertainties 
to the combined combination of Pq(ctv). 

For secluded DM models with a vector mediator, we showed in the previous section that 
the dominant primary product of annihilation will be e + e~ pairs, rather than 7-rays. While 
7's dont suffer the same energy degradation and so can have a larger observable flux due to 
annihilation throughout the halo, particularly when combined with secondary 7-production 
through bremsstrahlung and synchrotron emission, we will focus now on the observability of 
the positron component. As discussed earlier, the recent preliminary announcement from the 
PAMELA collaboration, apparently supporting the existing suggestive results from HEAT 
and AMS-01 of an excess positron fraction above 10 GeV, is difficult to reconcile with con- 
ventional WIMP models as it requires an annihilation cross-section well above that expected 
for a conventional thermal WIMP, i.e. (av) ~ 3 x 10~ 26 cm 3 s _1 . Thus, to explain the data 
an astrophysical 'boost factor' (p 2 ) / (p) 2 3> 1 is often postulated reflecting the possibility of 
a local enhancement in the number density due to substructure in the halo. 

Various analyses [HI E] indicate that the boost factor required to a produce a positron 




in 



dN e+ 
dE 




(20) 
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100 GeV 



1 TeV 



10 TeV 



100 GeV 



1 TeV 



10 TeV 



HI, 



Fi gure 2: The m^r^-a' parameter space for fermionic (left) and bosonic (right) secluded WIMPs, showing 
the cosmological abundance constraint (solid line), the region where radiative recombination into WIMP- 
onium is allowed (to the right of the dashed line), and a grey band where, limited by two extreme models of 
e + propagation, the positron excess would be above background for energies above 10 GeV. 

excess in the HEAT range, i.e. n e +/(n e - + n e +) ~ 0.1 above 10-20 GeV, is roughly: 



where N* is a multiplicity factor counting the number of positrons produced for each an- 
nihilation. In the present case, given the various decay channels through ortho- and para- 
states of WIMP-onium, we have iV* ~ 2. Our primary observation is that this boost-factor, 
although somewhat uncertain, appears to be of precisely the order of magnitude achieved 
through recombination of WIMPs in the galactic halo as discussed in the previous section. 
Thus, in such scenarios, one would expect an observationally detectable positron excess even 
in the absence of astrophysical boost factors. 

To illustrate the interesting region of parameter space, we plot in Fig. 2 the — 
a' plane taking the extreme value my ~ 4 MeV to maximize the allowed region where 
annihilation occurs via radiative recombination. This line intersects the constraint on the 
primordial abundance at around 200 GeV, and the region to the right is allowed. The grey 
band illustrates a representative range, allowing for unknown parameters in the propagation 
of positrons in the local region, that would produce an enhanced cross-section consistent 
with the observed excess seen by HEAT and AMS-1. The two extreme propagation models 
were taken from [27]. This band is almost parallel to the cosmological abundance line, but 
the enhancement factors J\f^^ lead to an overlap region which covers a large mass range 
due to the uncertainties in positron diffusion. Indeed, we see that there is broad consistency 
between the cosmological constraint and the apprent level of the positron excess for 200 GeV 
< < 10 TeV, and 200 GeV < < 4 TeV, given a value for m v consistent with (|2J). 
It is also worth noting that if my lies somewhat above the recombination range (j2J), the 
Sommerfeld enhancement factor ([Hj) still applies, and so the grey band can be continued to 
the left of the dashed line in Fig. 2, but shifted upward according to the reduction in the 
annihilation rate. 




where K 



0(100) 



(21) 
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4. Discussion 

In this note, we have presented a simple scenario within which the annihilation cross-section 
of a thermal WIMP may be significantly enhanced in the galactic halo with respect to its 
standard picobarn rate at freeze-out. For most WIMPs, the astrophysically relevant annihi- 
lation rate (av) is equal to or smaller than the corresponding rate at freeze-out. However, 
for some models this general rule is violated, and we discussed three generic mechanisms 
that could boost the halo annihilation signal: (A) (relatively) long-range interactions in the 
WIMP sector; (B) accidental near-threshold resonances; and (C) new annihilation channels. 
While (A) and (B) have been considered previously, we focused on the third possibility where 
the mechanism relies on the possibility of annihilation proceeding through recombination, via 
a metastable WIMP-onium bound state, and thus requires the presence of a light mediator 
and thus long-range forces in the dark sector. We observed that the recombinant rate is nat- 
urally enhanced by a factor of (9(100), while at the same time such models characteristically 
allow for a significant decoupling, or seclusion, of the dark sector with a small cross-section 
for scattering on nuclei. The vector nature of the mediator tends to enhance the leptonic 
decay fraction relative to 7's, which is particularly intriguing given recent speculation about 
a high-energy positron excess. While the excess observed by HEAT, AMS-01, and now in 
preliminary form by PAMELA, may find a conventional astrophysical explanation, it is in- 
teresting that a fairly generic class of WIMP models, namely secluded models with light 
mediators, can naturally have an enhanced annihilation cross-section in the halo that may 
have additonal consequences that would be interesting to explore further. We will conclude 
by mentioning several other aspects of secluded WIMP models. 

• Implications for BBN: It is well known that the annihilation of standard neutralino-like 
WIMPs during BBN does not affect the light elemental abundances unless the neu- 
tralino mass is under ~ 80 GeV [28]. However, the enhancement of the annihilation 
cross section resulting from a long-range interaction in the dark sector will change this 
conclusion as the number of WIMP annihilations after freeze-out acquires a different 
(milder) scaling with temperature. It is customary to characterize the impact of decay- 
ing or annihilating particles on the light elemental abundances by the quantity £ - the 
energy release in GeV normalized by the thermal photon number density. Assuming 
that all the energy from annihilation is released in the form of electromagnetic radia- 
tion, which is indeed the case for the models considered here, and following the earlier 
calculation in [23], we estimate that the value of £ due to annihilation at T ~ 0.1 keV 
can be as large as a fewxlO -13 , which will have an impact on the late-time generation 
of 6 Li at an observationally interesting level of 6 Li/H ~ 10~ n . 

• Signatures in the visible sector: Among the allowed 'portals', i.e. renormalizable cou- 
plings of the mediator to the SM: mixing of a U(l)' vector with the photon; coupling 
of a singlet scalar to the Higgs sector; or perhaps the right-handed neutrino, existing 
theoretical prejudice on the mass scale suggests that only the U(l)' can naturally be 
much lighter than the weak scale. When the scale of the mediator mass is much smaller 
than the weak scale but the coupling between the secluded and visible sectors is not 
too small, a new experimental signature may become available [31 [29]. For mixing 
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angles of order k ~ O(10~ 3 ) between the visible and secluded U(l)' sector, as would 
naturally arise from integrating out heavy states charged under both hypercharge and 
the secluded U(l)', the presence of 0(1-100 MeV)-scale mediators can create small 
but potentially detectable deviations from the SM in rare decays, the muon (g — 2) 
anomaly, and precison QED tests [30J. 

• Implications for direct detection and rare isotope searches: As noted in Sect. 2, the 
relatively long range force mediated by V bosons leads to a considerable enhancement 
of the direct detection rate, thus limiting the level to which the dark sector can be 
secluded. Its worth pointing out that if this effect is counterbalanced by the appro- 
priate splitting ArriM) of the WIMP states, it may lead to distinct signatures for the 
'inelastic' scattering of dark matter [T4]. At the same time, if m v l is much larger than 
the typical nuclear size of a few fm, sizable values for (e'tz) could possibly lead to the 
binding of WIMPs with heavy nuclei due to V^-exchange, resulting in recombination 
signatures in direct detection and the presence of stable anomalously heavy nuclei |16j . 
However, all inelastic scattering possibilities are quite sensitive to the parameters of 
the model. 

Note added: While this paper was being finalized, two preprints appeared [21] which also 
discuss the relevance of a long-range interaction in the dark matter sector, and overlap with 
some parts of the present paper. 
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